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The loop equations for a chain of hermitian random matrices are computed 



, expUcitely, including the 1/N corrections. To leading order, the master loop 

\ equation reduces to an algebraic equation, whose solution can be written in 

• terms of geometric properties of the underlying algebraic curve. In particular 

. we compute the free energy, the resolvents, the 2-loop functions and some 

mixed one loop functions. We also initiate the calculation of the expan- 
sion. 

^ '• 1 Introduction 



The multimatrix model, or chain of matrices is defined by the probability 
weight for J\f + 1 hermitean matrices of size N x N: 

Zj 

Z := (iV/27r)-^(-^+i)^' e-f^^ := / ciMo . . . ciM^ e'f E^-o v.(m,)-e1i a^.^a/,] 



(2) 

where (A; = 0, . . . , A/") are N x N hermitean matrices, and dM^ is the 
product of Lebesgue measures of all real components of M^. Z is the partition 
function, F is the free energy. T is the temperature, it can be chosen equal 



-'^Work partialy supported by the EC IHP Network: "Discrete geometries: from solid 
state physics to quantum gravity", HPRN-CT-1999-000161. 
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to 1 in most of the paper, except when one is interested in derivatives with 
respect to T. 

The potentials Vk are polynomials of degree dk + I: 



In order for Z to exist, we have to assume that all V^'s are bounded from 
below on the real axis. However, that constraint can be relaxed by studying 
"normal matrices" instead of hermitean matrices^, or by considering Z as a 
formal power series in the coefficients of the potentials. For simplicity, we 
assume here that the V^'s are real^. 

Each Mk lies in the potential well Vk and is linearly coupled to its neigh- 
boors Mk-i and Mk+i. Mq and Mj^ have only one neighboor, and thus the 
chain is open. So far, the problem of the closed chain has remained unsolved. 

The multimatrix model is a generalization of the 2-matrix model. It was 
often considered in the context of 2-dimensional quantum gravity and string 
theory. Its critical points are known to represent the minimal conformal field 
theories, characterized by a pair of integers {p,q)- It is known that one can 
get a {p, q) critical point with a multimatrix model where A/" = q — 2. The 
necessity of studying multimatrix models can be understood from the fact 
that the one matrix models contains only the critical models with q < 2. 

Recent progress have been made in the understanding of the two-matrix 
model [SlIUinilinillHllin], and it has been often noticed that the chain of 
matrices presents many similarities with the 2-matrix model |7j. 

The loop equations of the 2-matrix model have been known for a long time 
[T^l^ and have been concisely written in [TD^ including the 0(1/A^^) terms. 
The loop equations for the chain of matrices have been written in an appendix 
of and in a draft of J3] without the 0(l/iV^) terms, and without any 
proof, and also with misprints. The present paper is mainly a translation 
from french to english of the draft with proofs, and additional results. 
It thus fills a gap in the appendix of ^31 • In particular, the terms 
are included, opening the way to finding next to leading order corrections 

^i.e. the eigenvalues are located on complex paths instead of the real axis. For a 
potential of degree dfe + 1, there are dk homologicaly non equivalent complex paths on 
which Vk is bounded from below 

■^Most of the results presented here remain valid for complex V^s. The asumption of 
real potential merely allows to have a more concise derivation of loop equations. 




(3) 
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as in [in]. The leading order solution is also presented in an an updated 
algebraic geometry language. 

The paper is organized as follows: 

1. introduction 

2. definitions of the loop functions. 

3. the main result and its derivation: the master loop equation. Readers 
less interested in technical details can skip the derivation. 

4. to leading order, the master loop equation is an algebraic equation. We 
discuss its solution. 

5. computation of the free energy in the large N limit, along the lines 
drawn by 

6. computation of other resolvents and correlators in the large N limit. 

7. next to leading order, along the lines drawn by [TlU lllj. 

8. examples (2 and 3-matrix models, gaussian case, one-cut case). 

The most important results of this paper are the master loop equation 
eq.(I2Hl), the free energy eq. ()lll|) . the two-loop function eq. fjl62|) . the mixed 
correlators eq. ()165|) . eq. ()182p and eq. ()172|) . 

2 Definitions 

2.1 Moments 

We define the moments: 

Tno,n,,...,nj, ^ ( M^^^Mf^ . . . M^) + -L ( tr . . . M]^^M^'>) (4) 

where <> means averaging with the probability measure of eq.(0). If all 
the potentials V^'s are bounded from below, it is clear that the moments 
are well defined convergent integrals. If on the contrary the V^'s are not all 
bounded from below, then the moments are only formaly defined through 
their expansions in the coupling constants. 
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2.2 The complete one loop-function 



We define the complete one loop-function as follows: 

no,...,nj^=0 ^0 ^1 • • • ^AT 

which can also be written 

W{zo,...,zm) = (tr \- ^-ji---- ^-jr~/ 

2N \ zq- Mq zi - Ml zj^ - Mj^ I 

1/1 11 
H ( tr . . . 

2A^ \ ZM- Mm zi - Mizo- Mo 

Notice that it is only a formal definition (even when all the potentials are 
bounded from below), it makes sense only through its power expansion in 
the vicinity of all z^. oo. That function is merely a convenient concise 
notation for dealing with all the moments at once. 



2.3 Uncomplete one loop-functions 

We also define "uncomplete" loop-functions. For any integer k between zero 
and Af + I, and for any ordered subset of [0,Af] I = {ii, . . . , 4} {0 < ii < 
12 < . . . < ik ^ A/"), we define: 



^n,...,ifc (^ii 5 • • • 5 ^ifc) 2Af ( tr 



1 / +^ 1 1 1 



2N \ z,^-Mi^ ■ ■ ■ z,^-M,^ Zi^-Mi^ 

Res JJ dz.j W{zo, ...,zu) 
Pol W{z^,...,zm) T\z, (7) 
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where Pol means the polynomial part in the vicinity of oo. 



2.4 The resolvents 

As a special case of eq.(|7j), the resolvent of the fc*^ matrix of the chain is 
defined by: 

W^^(^^)-46r-^\ (8) 



N \ Zk-Mk 
4 



In particular, the resolvent of the first matrix of the chain will play a promi- 
nent role: 

The master loop equation derived in section EUl is an equation for Wo{zq) as 
a function of zq, and in section EJ we determine the large N limit of Wo{zq). 
In section lUTTl we derive the large N limits of the other resolvents W^- 

Instead of Wq{zq), it appears more convenient to consider the following 
function: 

Z,{zo) ■.= V^{zo)-TWo{zo) (10) 
and more generaly we define: 

Z_i(zo) ■.= TWo{zo) 

Zo{zo) ■■= Zq (11) 
Zk+i{zo) := Vl{Zk{zo)) - Zk-i{z^) l<k<M 

For short, we write: 

Zk := Zk{z^) (12) 

Notice that all Z^'s are polynomials in both zq and Wq{zq)^ and therefore, 
any polynomial of the Z^s is a polynomial in zq and Wq{zq). 



2.5 Polynomial loop functions P and U 

We define the functions fk,i{zk, . . . ,zi): 

fk,i := if / < A; - 1 

fk,k-l '■= 1 
fk,k{Zk) ■= Vl{Zk) 

fk,i+i{zk, . . . , zi+i) := Vl^-^{zi+i)fk,i{zk, ...,zi)- zizi+ifk,i-i{zk, • • • , ^/-i) 

(13) 

Then we define for k <l: 

Pk,i{zo, ...,zj^) := Pol fk,i{zk, • • • , zi)W{zo, . . . , z^f) (14) 

where Pol means the polynomial part in the vicinity of oo. Pk,i depends on 
A/" + 1 variables, and is polynomial in the variables Zk, ■ ■ ■ , zi. 
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By definition of fkj we liave: 

Pk,l{zo, ...,Zj^)= Pol Vi'{zi)Pk,l-iizo, ...,Zj^)- Pol Zi^i Zi Pk,l-2{Z0, ■■■,Zj^) 

(15) 

In particular, we define: 

P{zq, ...,z^):= Po,Ar(^o, ■■■,z^) (16) 

which is polynomial in all the variables, and 

U (^0, • • • , ^at) := Pi,Ar(^o, ■■■,z^) (17) 

which is polynomial in all the variables but Zq. 
We also define: 

E{zo, ...,zm):= (Vo'(2;o) - ^i)(V(r(^Ar) - ^AA-i) - TP{zo, ....z^) (18) 

which is also a polynomial in all variables. These polynomials will play an 
important role in the master loop equation. 

2.6 2-loop functions 

We define the loop insertion operator as in [T]: 

dVi{z)- zdgi,o ^^z'^+^dgi,k ^ ^ 

where Vi[z) = gi^ + J2k>i ^^'^i ^^cl the derivatives are taken at gi^k = for 
k > di + 1. It is a formal definition, which makes sense only order by order 
in its large z expansion. 

It is such that (for \z\ > \z'\): 

dVkjz') 1 dV^z') 1 

When applied to eq.(j2I), it produces a resolvent, i.e. the expectation of a 
trace: 

^ 4 (tr^) (21) 
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and more generaly, the action of g^jj^ on an expectation value inserts a new 



trace. 



In particular, we define the following two loop functions (i.e. two-traces): 



d 

W.i{zo, . . . , zm; z) ■= iy(2:o,...,2:^) 



1 



■tr 



2 \ zq- MqZi- Ml zj^- Mj^ z-Mi 

1/1 1 11 
+ - ( tr tr . . . 

2 \ z - Ml ZM - Mf4 zi - Ml Zq - Mo 

t2 



-N'Wi{z)W{zq,...,zm) 



(22) 



W, 



ii,...,ik 



;/(2^ii; • • • ; ^ifei z) . T 



d 



dViiz] 



"W^n,...,jfc(^ji; • • • ; ^ifcy 



\J\f+l-k 



Res Y\ dzj W.i{zo, ...,zj^; zf^S) 



HI 



Pk,i;jizo, . . . , zj^; z) := Pol fk,i{zk,...,zi)W.j{zo,...,zj^;z) (24) 

P-j{zo, ...,ZM\z) := Po,Ar;i(^o, ■ ■ ■ , zj^; z) (25) 
U.j{zo, ...,Zf^;z) := Pi^j^.j{zo, . . . , zj^; z) (26) 
In particular, the following two-loop functions will play an important role: 



Wk-i{z;z) = T- 



dViiz) dVkiz')dViiz) 

1 1 
tr — tr 



T 



dVJz' 



W,Jz; z') 



(27) 

z-Mi z'-Mk/ ^ ^ 

where the subscript c means connected part < AB >(..=< AB > — < A >< 
B >. 



3 Loop equations 



3.1 Main result: the master loop equation 

We prove below that the master loop equation (named after 
written: 



rp2 

E{Zo, Zi, . . . , Zjsf) = j^U-fl{ZQ, Zi, . . . , Zjsf] Zq) 



) can be 



(28) 



where the functions Zk{zQ) have been defined in eq. (jll|) . Examples for 
M = 1 (the 2-matrix model) and M = 2 (the 3-matrix model) are explicited 
in section IHl 

As an intermediate step in proving eq. ()28p , we need to prove the following 
formula for all 1 < < M: 



{Zk - Zk)WQ^k,...,J\f{zQ,Zk, . 


. . , Zf^) = Wo^k+l,...,Af{zo, Zk+1, . . 






■~-Po,fc-i(^0) Zi, . . . , Zk 


_1, Zfc, . . . , Zj^) 




~JpPl,k-l;o{Zo, Zi, . . . 


, Zk-1, Zk, ■ ■ ■ , zj^\ Zo) 



(29) 



3.2 Derivation of the master loop equation 

• Proof of eq.(j2ni) for A; = 1: 

The invariance of the matrix integral under the infinitesimal change of 
variable (see [T^ IHllTUj. and pay attention to the non-commutativity of 
matrices, and to the hermiticity of the change of variables): 

5Mo = . . . 

2 21 - Ml Z2 - Ms ZM- Mm zq - Mo 
T 1 1 11 

zo- MqZm- Mm'" Z2- Ma zi - Mi ^ ' 

implies 

T 

TWo{zq)W{zq, Zi,..., zm) + ^W^;o(2o, 2^1, • • • , ZM] Zq) 

= Vq{zo)W{zo, Zi,..., zm) - Po,o{,Zq, zi, . . . , z^) 

-ZiWqx,„^m{.Zq, Zi,..., Zm) + Wq^2,...,u{.Zq, Z2,..., Z^f) 

(31) 

The LHS comes from the Jacobian of the change of variable, and the 
RHS comes from the variation of the action. The first two terms of the 
RHS come from 5tr1/o(Mo) = tTV^{Mo)6Mo, and use of eq.®. The 
last two terms of the RHS come from 5tr MqMi = tr Mi^Mq, and use 
of eq.®. 

Using Zi = Vq{zo) — TWo{zq), this can be rewritten: 

{Zi - Zi)Wo^l^,„^MiZo, Zi, . . . , Zm) = Wo^2,...,Afi^O, Z2, . . . , Zm) 
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T 

~]^^;0*^^0' ^1' • • • ' ^at; Zo\?>2) 
Therefore we have derived eq.()29|) in the case k = 1. 

Proof of eq.(j2nD for A; = 2: 
Notice that eq. (jH^ imphes: 

- Res dzi W{zq, zi, 2:3, ... , zj^)Vl{zi) = V^i'(^i)W"o,2,...,a^(^o, 2^2, • • • , ^a^) 

"~-Po,i(^0) Zi, Z2, . . . , zj^) 

T 

— irp;Pl,l;Q{zo, Zi, Z2, . . . , Zjsf] Zq) 



(33) 



Ar2- 



The change of variable: 



2Z2-M2 z^-M^zq-Mo 2zo-Moz^-M^ Z2-M2 

(34) 

gives (since SMi is independant of Mi there is no jacobian, the LHS is 
zero, and the RHS is the variation of tr V^/(Mi) - MqMi - M1M2): 

= -Res dziW{zQ, zi, Z2, . . . , Zf^)V({zi) 

-ZoWo^2,...,Afi^O, Z2,..., Z^f) + W2,...,Afiz2, ■ ■ ■ , Z^f) 
-Z2Wo^2,...,Afi^O, ^2, • • • , Z^f) + VFo,3,...,A/-(^0, ^3, • • • , Zj^) 

(35) 

i.e., using eq.(|H^ and V({Zi) = zq + Z2: 

{Z2 - ^2)W^0,2,...,Ar(^0, ^2, • • • , Zj^) = lVo,3,...x(^0, ^3, • • • , Zj^) 

—Po,lizo, Zi, Z2, . . . , Zj\f) 

T 

~^^-Pl,1;o(^0) Zi, Z2, ■ ■ ■ , Zj^\ Zq) 

(36) 

Therefore we have derived eq.()29|) in the case k = 2. 
Proof of eq.(I2ni) by recursion on k: 
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We have already proved eg . (1291) for k = 1 and k = 2. Now assume 
that eq. (j29|) holds for k — 1 and k, we are going to prove it for + 1. 

Note that eg . (1291) for k — 1 implies (multiply by z^-i and take the 
residue at Zk-i oc and Zk — * oo): 

Res dzk-i Res dzk ^fc-iWo,fc-i,fc,..,Ar(^o, Zk-i, Zk, Zk+i, ...,zj^) 
= Zk-iWo^k+i,...,Af{zo, Zk+i, . . . , 

■^fc — 1 

T 



k 



(37) 

and eq.(j2ni) for k implies (multiply by Vl{zk) and take the residue at 

Zk oo): 

-Res dzk V^{zk) IVo, k,k+l,...,Af{Zo, Zk, Zk+l, . . . , 
= Vl{Zk)Wo^k+l,...,AfizO, Zk+1, ■.■,Zj^) 

- Pol Vl{zk)Po,k-i{zo, Zm)\,,=Zo,...,z,=Z, 
T 

Pol Vl.{zk)Pi,k^i-oizo, ...,zm; zo)\ 



(3J 



Then consider the change of variable: 
5Mk ^ ^ 



2 Zk+i - Mk+i zx - Mj^ zq - Mo 
111 1 



2zo~ MQZ_^f- Mj^ Zk+i - Mk+i 

it gives 



(39) 



= - Res dzk Vl,{zk) Wo^k,k+i,...,N{zo, Zk, Zk+i, . ■ ■ , zj^) 

- Res dzk-i Res dzk Zk-iWQ^k~i,...,Af{zo, Zk-i, ...,z^f) 

-Zk+lWQ^k+l,...,jv{zo, Zk+l, . . . , Zj^) + WQ^k+2,...,Af{zo, Zk+2, ■ ■ ■ , Zjsf) 

(40) 

using eq. ljHTj) and eq. (jHS|) . as well as VI,{Zk) = Zk-i + Zk+i, we get 
ea. d^ for k + 1. 
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eq.pUl for k = M: 

In particular, the previous recurrence derivation shows that eq. ()29|) for 
k = Af reads: 

{z^f - Z^f)Wo^j^{zo, zj^f) = Wo{zo) 

~ ]v2-Pl,A/'-l;o(^0) Zi, . . . , Zj^f_i, Zj^] Zq) 

(41) 

Proof of eq. (P|) : 

In particular, from eq. ()29|) for A; = A/" — 1 we derive: 

Res dzj^_iRes dz^f zm-iWo^^j^_^^zj^{zo, zj^-i, zj^) 
= Zm-xW^{z^) 

~ zM-hj^ ^^-i^Ar^o,Ar-2(2:o, • • • , ^Ar)L,=^„,...,,^=^^ 
T 

(42) 

and from eq. (|4H) . we derive: 

- Res dzj^ y'}^{zj^) VToxl^o, zm) 
= V^{Z^)Wo{zo) 

- Pol V^{zx)Po^x-i{zo, zx)\,,=Zo,...,z^=z^ 
T 

-Jp Pol VX-(^A')PlX-l;o(2:o, ■■■.ZM- ^o)L,=Zo,...,.^=Z^ 

(43) 

Then, the change of variables 



gives 



= - Res dzj^ Vj^{zj^) Wo,Ar(^o, ^at) 

- Res dzj^_i Res c^^at 2;Ar-iVFo,w_i,2Ar(^o, 2;Ar-i, ^at) (45) 

i.e., using eq. ljl^ and eq. (PH|l : 

= (VX-(^A^) - Zm-i)Wo{zo) 
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(46) 



Recalling that Z\ = Vq{zo) — TWq^zq), we get: 



(V^iZo) - Z,){V;^{Z^) - Z^_,) - TPox(^o, z,,..., Z^) 



-T7^Pl,Nfi{ZQ, Zi, . . . , Zj^] Zq) 



(47) 



which ends the proof of eq. p8|l . 

4 Large N leading order, algebraic equation 

Throughough this section and sections 5 and 6, we abusively denote with the 
same name, the loop functions and their large N limits. 

To leading order, the master loop equation reduces to an algebraic equa- 
tion for Wq as a function of Zq: 



where Zi = Vq{zo) - TWo{zo) and Zk+i = V^{Zk) - Zk-i- 

E is a. polynomial of given degrees in each variable, and with known 
leading coefficients. The problem is that most of the subleading coefficients 
of E{zq, zi, . . . , zj^) are not determined by the loop equations. They are 
determined by additional hypothesis, which will be explained in section 14.51 
Prior to that, we need to study the geometry of that algebraic equation. 

The complex curve Wq (equivalently Zi) as a function of zq is a one- 
dimensional submanifold of C-'^"'"^, which is the intersection of dimension 
M algebraic submanifolds. The curve Wq as a function of zq is thus a Riemann 
surface £. 

Instead of viewing the Z^s as (mutivalued) functions of zq, it is more 
appropriate to view the Z^'s as (monovalued) complex functions over £: 



E{zo,Zi,...,Zj^) = 



(48) 



(49) 



such that for all p E S: 



(50) 
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Since the function Zo{p) is not injective, i.e. the point p such that Zo{p) = zq 
is not unique, the Z^s are multivalued functions of zq. Similarly, one could 
consider any Zk as a multivalued function of any zi. 

Instead of dealing with multivalued functions, we slice £ into domains 
called sheets, such that in each domain the functions we are considering are 
bijections. The 2;o-sheets are thus domains in which the function Zq{p) is 
injective. More generaly, the 2;fc-sheets are domains in which the function 
Zkij)) is injective. 

Let us study the zo-sheets first. For that purpose, we define: 



(51) 





= -1 




= 1 


Tk 


= dodi . 




:= -1 


SM 


:= 1 


Sk 


:= dk+i 



4.1 ZQ-sheets 

The equation 

Eizo,Z,,...,Z^) = (52) 

has degree 1 + did2 . . . dj^ in Wq and degree d^di . . . dj^ in zq. Therefore, 
W^o(-2^o) (or equivalently Zi) is a multivalued function of zq which takes 1 + 
did2 . . . df^ values: in other words, there are 1 + did2 . . .d^f ^o-sheets. We 
identify these sheets by looking at the asymptotics of Zj^ when zq ^ oo. 

• The physical sheet 

From eq.©, there must exist at least one solution of the algebraic 
equation such that: 



which implies that: 
and in particular 



Woizo) + 0{l/zl) 

20-*oo Zq 



Zk 



(53) 



(54) 



z;^ = 0(;2o"°'^^-"^-^) (55) 

The Zo-sheet in which these asymptotics hold is called the 2;o-physical 
sheet. 



dodi...d^_ 
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• Other sheets 

Notice that the equation for Wj\f{z^f) as a function of zj\f is the same 
algebraic equation. In other words, there is a solution of the algebraic 
equation which is such that: 

Zk = 0{zp) , Sk = 4+1 • • • dx_idj^ (56) 

= 0(^;/'^-'-^-'^) , z, = 0{zy (57) 

Since the number of s^q roots of unity is exactly Sq = di . . . dj^_idj^, 
we have all the solutions. 

We thus have rQ + sq = 1 + di . . . dj^^idj^ sheets. In one of them we have 
ZM ~ 0(zq-^), and in the di. . . d^^idj^ others we have zq ~ 0{z^). 

4.2 Algebraic geometry 

We now consider the algebraic curve £^ in a more geometric language. An 
abstract point p E £ can be represented as a couple {zq, Zi) such that zi = 
Zi{zo), or by any other parametrization. For instance it can be described as 
a point of C-'^, at the intersection of Af codimension 1 manifolds. 

Algebraic geometry is an active and important part of mathematics, and 
lots of tools have been invented to describe the geometry of algebraic Rie- 
mann surfaces. We refer the reader to fOJEI for an introduction. 

A Riemann surface is locally homomorphic to the complex plane C, that 
means that small domains of S can be maped on small domains of C. The 
map, which is one-to-one and holomorphic in that domain is called a local 
parameter, let us call it: 

pes ^ x{p) e C (58) 

Any complex valued analytic function on S can be localy represented by an 
analytic function of x{p). 

For instance, Zq = zq{p) is often a good local parameter on £, except 
when Zq approaches a singularity. Therefore all functions on S can be localy 
written as functions of zq. They can also be locally written as functions of 
any Zk. 

Notice that the function Zo{p) is not injective, it takes the same value 
Zq for different points p, namely the same value of Zq corresponds to 1 + Sq 
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points p. Therefore the function Zi{p) which is a well defined monovalued 
function on £^ is a multivalued function of Zq. The Zo-sheets are domains of 
S where the function zo{p) is injective. In particular in each sheet, there is 
only one point Poo where zq ^ oo. 

Let poo+ be the point in the physical sheet, such that zo{poo+) = 00. In 
the vicinity of Poo+, we have eq.(|^: 

Zk{p) ~ ZQ^ip) (59) 

i.e. Z}.{zq) is analytical, which indicates that Zq is a good local parameter 
near p^+. 

Now, let Poo- be the point at 00 in the ^A^-physical sheet, i.e. where 
zm{Poo~) = cxD and with behavior eq. (j^ : 

Zk{p) ~ ^^(p) (60) 

i.e. Zk{zj\f) is analytical, which indicates that zj\j- is a good local parameter 
near Poo~- In the vicinity of Poo-, the function Zq{p) takes the same value zq, 
Sq times, therefore there are Sq sheets which meet at Poo-- It is clear that zq 
is not a local parameter near Poo-, but zj^ is. 

It is also clear that there can be no other point p such that zq{p) = 00. 
The algebraic curve S has only two points at 00. 

Notice that the intermediate z^^^s with Q < k < M are not appropriate 
local coordinates near Poo+ neither near Poo- (unless many of the rf^'s are 
equal to 1). 

Let us summarize as follows: 

- The function Zk{p) has a pole of degree = d^di . . . dk-i near Poo+- ^o(p) 
is a local parameter near Poo+- 

- The function z^ip) has a pole of degree = d^+i ■ ■ ■ dj^^idj^f near Poo-- 
zj^{p) is a local parameter near Poo-- 

- The function Zk{p) has no other pole. 

Remark: From eq. ljSHj) . it is easy to prove by recursion that: 
Res Zk-idzk = — T = — Res Zk+idzk , Res Zk+idzk = —T = — Res Zk-idzk 

P00+ P00+ Poo- Poo- 

(61) 
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4.2.1 Genus and cycles 

Let g be the genus oi £. It can be proved (using the Riemann-Roch theorem, 
or the method of [12]) that 

M 

g < := J]^ 4 - 1 (62) 
fe=o 

Notice that the polynomial P has 1 + (^max coefficients, and its leading coef- 
ficient is fixed, i.e. gra&x is the number of coefficients of P not fixed by the 
loop equations. 

Let Ai, Bi{i = 1, . . . , g) he a canonical basis of non-trivial cycles on S: 

Aif]Bj=5ij (63) 

the choice of non-trivial cycles is not unique, and we will see below one 
possible convenient choice. 

4.2.2 Endpoints 

If Zi is considered as a function of z^. {k ^ I), it has singularities (branch 
points) everytime that: 

dzkip) = (64) 

Indeed at such a point we have: z[{zk) = dzi/dzk oo. 

The zeroes of dzk{p) are called the endpoints, and are noted: 

Cfc^j , (i = 1, . . . , rfc + Sfe + 251) (65) 

Generically, the zeroes of dzk are simple zeroes and they are all distinct, 
which means that Zk behaves as a quadratic function of a local parameter 
X, while zi behaves linearly in x. Therefore, zi{zk) has a square root branch 
point near Zk{ek,i). 

4.2.3 Critical points 

It may happen that some endpoints coincide. This is called a critical point. It 
is not a generic situation, it happens only if the potentials V^'s are fine tuned 
to some critical potentials. The critical points are relevant for finding the 
representations of {p, q) minimal conformal models. A (p, q) critical point, is 
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such that there exist k ^ I, and a point e & £ such that: dz^ has a zero of 
degree p — 1 at e and dzi has a zero of degree g — 1 at e. Near e, 2;^ as a 
function of zi behaves with a,p/q exponent. 

From now on, we assume that the potentials Vk are not critical, i.e. the 
endpoints are all simple and distinct. 



4.2.4 Cuts 

The cuts in £ are the contours which border the sheets. Viewed in the 
Zfc-plane they are lines in C joining two endpoints. 

That choice of lines (i.e. the choice of sheets) is largely arbitrary, provided 
that they indeed go through the endpoints. 

There is a canonical choice, which comes from [13J, given by the contours 
where the asymptotics of the biorthogonal polynomials and their Fourier 
and Hilbert transforms, have discontinuities, i.e. some Stokes lines. That 
canonical choice is defined as follows: the set of Zk-cnt?, (contours on £) is 
the set of points p & £ such that there exists p' ^ p with Zk{p') = Zk{p) and 



3ft / Zk+idzk = (66) 




remark that for all p and p' such that Zk{p) = Zk{p') we have: 

p' rp' rp' 

Zk+idzk = / (Vk^Zk) - Zk-i)dzk = - Zk-idzk (67) 
J p J p 

so that the canonical cuts are left unchanged if we reverse the order of the 
chain. 



4.3 Sheet geometry 

We are going to describe briefly the sheet geometry, see fig. ()4.3p for a better 
vizualization. 

The equation Zk{p) = Zk has + Sk solutions for p, i.e. the curve £ is 
divided into + Sk sheets. When Zk ^ 00, Vk solutions approach Poo+, while 
Sk approach Poo-- In other words, sheets contain Poo+, while Sk sheets 
contain Poo-- 

Let us denote Ck the contour which separates the reunion of sheets con- 
taining P00+ from the reunion of sheets containing Poo-, oriented such that 
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it turns around Poo+ in the positive direction. That contour Ck will play an 
important role later. Note that is not necessarily connected (as in our 
example fig. 14. 3p . 
Let us also denote: 

• p+j,fc(zfc) (j = 1, . . . , Tfc) all the solutions of Zk{p) = Zk which are on the 
same side of than Poo+5 

• p-j^k{zk) (j = 1; • • • ; Sk) all the solutions of Zk{p) = z^ which are on the 
same side of Ck than Poo-- 



4.4 Filling fractions 

It is well known that the discontinuities of Wq{zq), i.e. the discontinuities 
of Z_i{zq) along the Zg-cuts in the physical sheet are related to the large N 
average density of eigenvalues of the matrix Mq: 

Poizo) = -7^{Wo{zo + tO) - Woizo - ^0)) (68) 

The support of po is the set of cuts of Z_i(zo), with endpoints the 2;o(eo,j)'s 
which belong to the physical sheet. 

The ratio of the average number of eigenvalues in a given connected com- 
ponent of the support of po to the total number of eigenvalues is a contour 
integral: 

T^^i = / p{zo)dzo = -77-^ i z_idzQ (69) 



T Jzo{eo,i) 2i7rT „ 

where the contour of integration on £ is one of the 2;o^cuts defined in section 
14.2.41 oriented in the clockwise direction, is non-zero only if the contour is 
a non-trivial contour, or if it encloses a pole Poo+ or Poo-- 

Therefore, there are at most g + 1 possible values of e^, the support of the 
density has at most g + 1 connected components. 

It is possible to choose the potentials Vi, . . . , Vj\f such that all the non 
trivial cycles Ai are in the physical sheet, and such that the cycles Ai are all 
cuts. The filling fractions in that case are the ^-cycles integrals: 

z_idzQ (70) 

For more generic potentials, we define the filling fractions by eq. ()70|) . 
eventhough they don't really correspond to numbers of eigenvalues. The 
cuts are integer homological linear combination of ^-cycles. 



2m 
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Figure 1: Sheets for Zk. example with = 6, = 4 and g — 1- The two 
black cirlces represent P(X)+ a-nd Poo-- The empty dries are the Tk + + 
2g endpoints. The thick lines, going through the endpoints, are the cuts 
separating the sheets. Each region separated by lines is a sheet, there are 
f^k + "Sfc sheets. The sheets form "flowers" near Pao+ and Poo- (the petals 
have an angle 27r/rfc near poo+ and 2ir/sk near Poo-)- There are g "handles", 
with A and B cycles: here A can be chosen as a cut, and B is represented 
as a thick dash-dot line. Each cross x (one in each sheet) represents a point 
p & £, such that Zkip) — z^. 
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We also define the chemical potentials as the S-cycle integrals: 



r. :=- 




Z-idzi, 



■0 



(71) 



We will see in part 16. II that this description remains valid for the densities 
of all matrices Mk (and not only Mq). We will prove (see eq. ljlfiOj) ) that Pk{zk) 
the density of eigenvalues of Mk is the discontinuity of Zk+i (or equivalently 
Zk-i) along the cuts. 

It is easy to prove, by integration by parts, that: 



If we anticipate on section 16.11 this proves that the filling fractions are the 
same for all densities, up to integer linear combinations. If the support of the 
density of eigenvalues of has < g + 1 connected components [a^^j, bk,i] 
{i = 1, . . . , mk), the filling fractions in each connected component are: 



where the coefficients Ak^ij are integers (possibly nul or negative). 
4.5 A or B cycles fixed ? 

As we said before, the loop equation is not sufficient to determine all the 
unknown coefficients of E. Some additional hypothesis are needed, two of 
them are often considered in the litterature: 

• condition B: if all the potentials are bounded from below, and the 
partition function is well defined, one is interested in finding large N 
limits for the free energy and various expectation values of traces of 
powers of the matrices. In that approach, one has to find a solution 
of the loop equation which gives an absolute minimum of the free en- 
ergy. The genus g and the filling fractions ej are not known, they are 




(72) 




(73) 




(74) 
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determined by the minimization condition dF/dti = 0, which reduces 
to: 

= 1, . . . , g (p zidzQ = (75) 

All the B cycle integrals must vanish. That condition is sufficient to 
determine g and all the e^, and it is sufficient to determine the unknown 
coefficients of E. 

This allows to find the large N limit of the free energy. Subleading 
large N corrections, i.e. the so-called topological large expansion 
of the free energy, exist only if g = 0, as was shown in jHlCD]. If the 
genus g is > 1, there exist asymptotics for the subleading corrections 
to the free energy, which have no 1/iV expansion, they are oscillating 
functions. 

• condition A: if one is only interested in the formal large A^ expansion 
of the free energy, one has to find a solution of the loop equation which 
corresponds to the perturbation of a given local extremum of the free 
energy. Therefore, g and the e^'s are fixed parameters which character- 
ize the minimum around which we perform the perturbative expansion 
(moduli). The following equations 

Vi = 1, . . . , (7 ® Zidzo = 2iT[ei (76) 

are sufficient to determine all the unknown coefficients of E. 

In the following, we will always assume that we have condition A, unless 
specified. Most of the results obtained for condition A, immediately translate 
to condition B, by exchanging the roles of A and i3-cycles, like the large A^ 
free energy. But some quantities, like the large A^ 2-loop functions, get 
contributions from the oscillating asymptotics, and will not be computed in 
this article. So, from now on, we assume condition A. 

4.6 Parametrization 

It is possible to parametrize the algebraic curve S in terms of 6 functions. 

Given g and a canonical basis of non-trivial cycles Ai, Bi{i = 1, g) , it 
is known [inillTj that there exists a unique basis of holomorphic 1-forms dui 
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on 8 such that: 

(h duj = 5ij (77) 
The matrix of B periods is defined as: 



Tij := (h duj = Tji (78) 
Jb, 

Tij IS cL symmetric g x g matrix, with positive imaginary part. 

If we choose an arbitrary basepoint po G S, we define the Abel map: 

rp 

u{p) := / du (79) 

which defines an embedding of S into C^. 

For any < A; < M, Zk{p) is a function on £, with a pole of degree at 
Poo+5 and a pole of degree Sk at Poo-- Zk must have + Sk zeroes, which we 
denote Pfc,j(0), (z = 1, . . . , + s^). The zeroes must satisfy: 

M(Pfc,i(0)) = rku{poo+) + Sku{poo-) (80) 



1=1 



then: 



^,(p) = YyiV'mp)-u{Pk,m)-z-r) 

^ ^6{u{p) - u{poo+) - z\ tY'=6{u{p) - u{poo-) - z] tY^ 

where 9 is the Riemann theta function [TBlll7j . z is an arbitrary non-singular 
zero of 6. The ratio of the RHS and the LHS is a function on £ with no pole, 
therefore it must be a constant, which we call Ak. 

We have a parametrization of £ in terms of 9 functions. 



4.7 Two-loop functions and Bergmann kernel 
4.7.1 The Bergmann kernel 

The Bergmann kernel B{j)\ p') is the unique bilinear differential form on £^£ 
with the following properties: 
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• B{p; p') , as a function of p, is a meromorphic form, with only one double 
pole aX p = p' with no residue, and such that in any local parameter 
z{p) we have: 

5. n ^ dz{p)dz{p') 

•y^ = l,...,9 §^^^^B{p,p')=0. 
4.7.2 The 2-point function 

For < k < J\f and < / < Af, define the following meromorphic differential 
forms on £ X C: 



dzk-i{p) 



dVi{z) 
dzk+i{p) 



dVi{z) 



dzk ip) dz 



dzk ip) dz 



We have the "thermodynamic identity" : 

Bk,i{P'^ z) = -Bk+i,i{p; z) 
and Zk-i + Zk+i = V^{zk) implies: 

dzi (p) dz 



Bk,iip; z) + Bk,iip; z) 



z-zi{p)y 



(83) 
(84) 

(85) 
(86) 



We are going to prove below that: 

BkAp-, z) = - E;Li B{p;p.,j{z)) 0<k<l 

Bk,i{p-, z) = E;Li B{p; p+,,iiz)) l<k<U 

BkAp-^ z) = E;Li B{p;p^,Az)) 0<k<l 

BkAp-^ z) = - E?=i B{p;p+,A^)) l<k<U 



(87) 



where B is the Bergmann kernel, and the Pj^'s have been defined in 14.31 



Proof of eq.dHZD 

First notice that eq. (jH3|) and eq. ljHUI) imply: 

BkAp; z) = BoAp; z) o<k<i 
BkAp; z) = -Bj^Av^ l<k<M 



(88) 
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and: 

[z - Zi[p)Y 

Consider the two-loop functions introduced in eq.(j2Zj): 

VTo z) = { tr — tr — ) = — ^- ' ,, , (90) 

W_\f-i(z; z ) = { tr — — tr — ) = - — , ,, , , (91) 

' \ z-M^ z'-Mi/^^^^ dzM{p-iA^))dz> ^ ' 



The expectation values tr jzjj^ tr y^j^^j and tr j^j^ tr pzjg^y are well 
defined integrals for ^ in the physical sheet, therefore Bq^i{p; z') can have 
no pole when z' = zi{p) if p is in the Zo-physical sheet, and Bj^^p^z') can 
have no pole when z' = zi{p) if p is in the 2;_/\^-physical sheet. Note that 
the derivatives d/dVi{z') are formally defined only for large z', which implies 
that: 

Bo^i{p; z') can have no pole when p = p+j^i{z'), j = 1, . . . ,ri, 
Bj\f^i{p; z') can have no pole when p = p^j^z'), j = 1, . . . , S/. 

Since the RHS of eq. (jHIJ|) has poles at all p = p^j^z'), we must have: 
Bo,i{p; z') can have no pole when p = p^j^z'), j = 1, . . . ,ri, 
Bq^i{p; z') has double poles when p = p^j^z'), j = 1, . . . , s;. 

That implies that: 

Sl 

Bi{p-z') ■.= Bo,i{p;z') + Y,Bip;P-j,i{^')) (92) 

i=i 

has no pole when zi{p) = z'. 

It obeys the following properties: 

• Since PFo(^o) behaves as I/zq near oo in the physical sheet, Wq-i{zo; z') 
must behave as 0{l/z'^) near Poo+, while dzo has a double pole at Poo+- 
Therefore -Bo,«(P! ^0 has no pole when p Poo+, and thus Bi{p; z') has 
no pole when p ^ Pc 



Similarly, Bj^^p; z') has no pole when p —>■ Poc^, and with eq. ljH^ . that 
implies that Bq i{p; z') has no pole when p Poo-, and thus Bi{p; z') 
has no pole when p — > Pc 



24 



• Near an endpoint e^^j, Zk-i{zk) has a square root singularity, i.e. the 
derivative dzk_i{p) /dVi{z')\z,,{p) has an inverse square root singularity, 
i.e. a simple pole which is exactly compensated by the zero of dzk{p)- 
Therefore Bk,i{p; z') has no pole when p — > ek,i, and so for Bi{p; z'). 

• Near any other point, Zk-i is an analytical function of Zk-, and thus Bk,i 
is analytical. 

• Bk,i must satisfy: 

i Bk,i{p;z') = (93) 

and so must Bi. 

Finaly we find that Bi is a meromorphic form with no pole, with all its 
A cycle integrals vanishing, therefore: Bi = 0, QED. 

Remark: 

We have: 

B{p, p') = dpdp, In e{u{p) - u{p') - z) (94) 

Notice that 

du,{p') = ^l B{p,p') = dz,{p')^ J^] (95) 
^^'^ JpeB, 2m dVo{zo{p'}) 

Once B is known, eq. ()95|1 give an explicit way of computing the Abel map 
u{p). 

4.8 Abelian differential of the 3rd kind and tempera- 
ture 

4.8.1 Abelian differential of the 3rd kind 

There is a unique abelian differential dS on S with only two simple poles at 
Poo±, such that: 

Res dS = - Res dS = 1 , Wi = l,...,g (p dS = (96) 

POC.+ Poo- /_4^ 

It has the property that: 

dS{p)= / B{p,p') (97) 
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We define: 

r]i := / dS (98) 

And, given a basepoint po, we define tlie following multivalued function 
on E: 

A{p) :^e-^po'^^ (99) 

A has a simple pole at Poo+ and a single zero at Poo-, therefore the following 
quantities are well defined: 

7 := lim ^ , 7 := lim z^{p)A{p) (100) 

Notice that the product 77 is independent of the choice of the basepointpo- 

4.8.2 Derivatives with respect to T 

Consider the abehan differentials: 

dzk+iip) 



dzk{p) , dSk{p) :-- 



Zkip) ^ 



dzk{p) (101) 



dT 

We clearly have (from Zk+i + Zk-i — V^{zk) and thermodynamic identity): 

dSkip) = -dSkip) = dSk+iip) (102) 

Therefore dSk and dSk are independent of k. 
They have the following properties: 

• near Poo+, -Z-i ~ T/zq + 0{zq'^), thus dSo ~ dzo/zo, i.e. dSo has a 
single pole at Poo+, with residue —1. 

• Similarly liSV has a single pole at Poo-, with residue —1, i.e. dSo has 
a single pole at Poo-, with residue +1. 

• near an endpoint ek,i, Zk-i{zk) has a square root singularity, thus 

has an inverse square root singularity, i.e. a single pole at 

Cfe j, which is compensated by the single zero of dzx, and therefore dSo 
has no pole at ek,i- 

• near any other point, Zk is a local parameter and Zk+i is analytical in 
Zk, i.e. dSo has no pole. 
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• The Ai cycle integrals of dSo vanish. 

There is a unique abelian differential with such properties, it is the abelian 
differential of the third kind defined in eq. (|96j) : 

dSkip) = -dSkip) = dS{p) (103) 

In particular, we have that: 

^^^ = V^, ^e, = (104) 

4.9 Derivatives with respect to Ci (condition A) 

Consider: 



duk,i[p) ■■= — 



2m den 



dzk[p) , duk,i[p) ■■= — 



dzkip) 

Zk{p) 

(105) 

We clearly have (from Zk+i + Zk-i = V^{zk) and thermodynamic identity): 

duk,i{p) = -duk,i{p) = duk+i,i{p) (106) 

Therefore duk^i and duk,i are independent of k. 

Following the same lines as in the previous section, we show that duk^i 
and duk^i have no poles, i.e. they are holomorphic one-forms. Moreover we 
have: 

de- 

duk,^ =0^ = kj (107) 

There exists a unique set of holomorphic one-forms with those properties, it 
is the holomorphic forms dui introduced in eq. ()77|) . therefore: 

duk,i = -duk,i = dui (108) 

This shows that: 

—^ = 2t7Tnj 109 
dei 

where Tij is the matrix of periods introduced in eq.((ZHl)- 
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5 Large free energy 

5.1 The large free energy 

The free energy F defined in eq.© has a large limit: 

F = + 0{N-^) 
We prove below, a generalization of the formula of [3]: 



(110) 













-T2(l+ArinT) 






= Res {Vk{zk) 

^ ' Poo 

fc=0 


- ^ZkVl.{zk))zk-idzk 




-T^{l+M\nT) 



;iii) 



where /i is the generalized version of [3^ and is defined as follows for any 
/i := / ( z^i)dzo+ / ( zj^+i)dzj^ -T\nzo{p) -T\nzj^{p) 

Jpoo+ ^0 Jp^_ ZM 

N Af 

+^yk{zk{p)) -yzk-i{p)zk{p) 
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k=0 



k=l 



H is independent of p (indeed d^ is a telescopic sum which cancels com- 
pletely). 

Proof of eq. (jllip : We remind that we assume condition A. We define: 

^ 1 
AK := 2T/i + 2 V qF, - 4F(°) + V Res {Vk{zk) - -ZkZk+i)zk+idzk 

+ Res {Vk{zk) - -ZkZk^i)zk-idzk (113) 



fc=0 



Notice that ReSp^_^ ZkZk-iZk+idzk + ReSp^_ ZkZk-iZk+idzk = 0, so that the 
expression in eq. ()113j) is the same as in eq. ()lll|l . 
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Let us compute: 
dK 



dVi{z) 



+ y^Res 



dVi{z] 
dzk+i 



+ ^Res 



poo+ dVi{z) 
dzk-i 



dVi{z) 



iVk{zk) - ZkZk+ijdzk 



+2T- 



Poo- dVi{z) 



2fc 



(Vfc(2;fc) - ZkZk-ijdzk 



we introduce the following multivalued functions: 



fc+i 



dVi{z) 



dzu 



Ckip) :-- 



2fe 



dVliz) 



dzk 



Zfc 



By integration by parts, and using eg. ()2()j) . we have: 

4-^- = - Res -^^dzi - Res -^^dzi + ATWAz) 

oVi[z) Poo+ Z — Zi Poo- Z — -2/ 



Res Ck{{Vl{zk) - 2fc+i)c^2fc - Zkdzk+i) 

, Poo + 

ERes Cfc((K(2;fe) - Zk-i)dzk - Zkdzk-i) 
n 



+2T 



9yU 



Res - Res -^^dzi + 4TWi{z) 

Poo+ Z — Zi Poo- Z — 2; 

Res Ck{zk-idzk - Zkdzk+i) 

, Poo + 

Res Ckizk+idzk - Zkdzk-i) 

Poc — 

k 



+2T +2y 

— Res 



-(i^;; — Res 

Poo+ Z — Zi Poo- Z — -Z; 



dzi + 4:TWi{z) 



(114) 



(115) 



+ Res Cj^zj^dzj^+i - C,-iz_idzo - V'(Cfc - Ck-\)zk-idzk 

n I ' ^ 



k=0 
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+ Res C,oZodz_i - C,M+\ZM+idz^ - y'(Cfe - C,k+i)zk+idzk 

Poo— 



+2T 



dfj, 



dVi{z) 



'dVi{z) 



k=0 



(116) 



Near poo+ we have z-i ^ T/zq + 0{zq therefore zq ~ T/z-i + 0(1), and 
thus C-i has a zero at Poo+- That implies that: 



Res (-iZ-idzQ = , Res (j^^iZj^/^idzf^ = 

Poo+ Poo- 



(117) 



Moreover, notice that 

Cfe(p)-Cfe-i(p) = 



(^) 



c^^;/;-! — — 



dzk- 



zk-1 



dVliz) 



dzk, therefore: 



Poo-\ 



dVi{z) 



+ 



dzk-i 



dVi{z) 



dzk 



Zk, 



dViizk] 



dzk 



kd 



dVi{z) zi{p)-z 
(118) 



Thus: 
, dK 



dVi{z) 



Res ^J=^^l±^dzi + Res ^J±^^^dzi + ATWi{z) 

Poo+ Z — Zl Poo- Z — Zl 



dfj, 



+ Res CAf^ArdzAr+i + Res Co^^oc^^-i + 27"- 

Poo+ Poo- OVl{Z) 

Res ^hl^^dzi + Res ^'^' ~^'-' dzi + 4:TWi{z) 



Poo4 



z- Zl 



Poo- 



2; - 2;; 



djjj 

Res CAr^AAC?^A/-+i - Res + 2T 

Poo- Poo+ 0*1// (2;) 



^V^(z) 



~ XI 2^ j[ ^^^^Si (Co^oo?^-i) - XI 2^ j[ ^^^^S' (C^/^ATfi^AA+i) 



:il9) 



where we have used Riemann bihnear identity, and disc means taking dis- 
continuity of the considered multivalued function when crossing cycles. 

Now compute d^/dVi{z), all terms cancel but: 



dVi{z) z-zi(p) 



dz_i 



dVi{z) 
30 



dzn — 



dVi{z) 



dzu (120) 



which is independent of p. In particular for p = Poo±, this proves that Co(Poo+) 
and Ca^(Poo-) are finite, and: 



dVi{z) 



Co(Poo+) = C^iPoo-) 



;i2i) 



Thus: 

Res Cozodz-i = TC,o{po 



djjj 

T r,,,, ^ = TCAr(Poo-) = Res CAfZ^fdz^f+i 
oVi(z) Poo- 



(122) 

Now compute the discontinuities (crossing Bi is equivalent to going around 



A): 



discB^ Co 



dz^i 



dVi{z) 



dzo = —2iTT 



d€i 



dVi{z) 



this implies that Co has no discontinuity along Sj, and 



disc^, Co 



dz_i 



dZn 



dVi{z) 



(123) 



(124) 



this implies that Co has a constant discontinuity along Ai. Moreover zq and 
dz-i are monovalued, i.e. they have no discontinuity along Ai or Bi. Thus: 



disc^^ {Cozodz^i 



dVi 



A, 



Zodz^i = —2mei 



A, 



dVi{z) 



Using similar arugments for zj^, we arrive at: 



dK Zl_l-Zl + l Zl + l-Zl^l .rr^.rf ^ 

4- — — = Res —dzi + Res — dzi + ATWi{z) 

oVi{z) Poo+ Z — Zi Poo- z — Zi 

When / = 0, by definition eq.(fTT|). we have 



(125) 



(126) 



TWq{z) = Res zidzo = — Res Z-idzo 

Poo+ Z — Zq Poo+ Z — Zq 

= Res z^idzo = — Res Zidzo 



therefore: 



Poo- Z — Zq 

dK 



Poo- Z — Zq 



dV^iz) 
31 



(127) 



;i28) 



which proves that K is independent of Vq. In particular we can choose Vq 
quadratic, and then we integrate Mq out, i.e. we reduce the problem of a 
chain of lenght M with potentials Vq, . . . , to a chain of lenght M —1 with 
potentials Vi — • • • > ^M- K is easy to check directly that: 

FmC^, V,,...,Vm) = - i, . . . , V^) - ^ InT (129) 

and thus: 

Ku{V,,...,VM) = KM{^,V,,...,VM)=KM-i{V^-^^,...,VM) + ^\nT 

(130) 

Therefore, by recursion on A/", we find that K is independent of all Vs. K 
could still depend on T and the ej's. K has been computed for A/" = 1 and 
A^ = [3", and we find: 

K = ^{l+U\nT) (131) 
which is the same result as if all potentials are chosen gaussian (see section 
QED 

Remark: En route, we have proved that for all /, the resolvent of the 
matrix is: 

Wi{z) = ^Res — - — zi^idzi = -^Res — - — Zi^idzi 

I P00+ Z — Zi I P00+ Z — Zi 

= —Res zi-idzi = -—Res zi+idzi (132) 

I Poo- Z — Zl I Poo- Z — Zl 



Remark: we have: 

/AT N \ 

/i = -Tln77 + ^e^r^i + Res \'^Vk{zk) -^Zk-iZk\ dS (133) 

i \fe=0 k=l / 

where dS and 7 and 7 and rji are defined in eq.ljHEl) and eq. ()l()()|) . 
Proof of eq. (|133p : Consider the functions: 

Mp) ■= ( z^i)dzo , (t)jv{p) := / ( zj^+i)dzj^ (134) 

Jpoo+ ^0 Jp^_ z^f 
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They satisfy: 

Res (podS = , Res (pj^dS = 
Poo+ Poo- 
following discontinuities along Ai or Bii 

= Ti , discg^ 00 = -discg^ (pj^ = 2m e 

have: 



(135) 



and they have the 

disc^^ 00 = -disc^. 
From eq. ()112|) and eq.(j96J, we 

= Res ^dS = Res (podS + Res (pj\fdS — T In 77 

0+ Poc+ Pao + 

/AT Af \ 

Res I ^ Vk{zk{p)) - ^ Zk-i{p)zk{p) j dS 

\fc=0 fc=l / 

use Riemann's bilinear identity 



(136) 



+ : 



(137) 



We need to compute ReSp^^ (pjvdS, we 

Res 4>^fdS = — Res (pj^dS + ® disc^- (pj^fdS — ^— ® discg. (pj^dS 

P00+ Poo- ^-^ 2m J A. ^-^2mjK.. 



= XI ^^^^ 



(138) 



QED. 



5.2 Derivatives with respect to 

We have: 



indeed, from eq. pi2|l . we have 
9ei 



+1 



Poo- 



dzj^ 



(139) 



(140) 



using eq . (11081) we have: 



^ = 2m I dui = r]i 



(141) 



We have: 



(142) 
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Indeed, from eq. (|lll|) . and using eg . (11081) . we have: 

4— — = 2Trii + 21^ + 2 > ejn^ + > Res {Vk{zk) - ZkZk+i)dui 

ERes {Vk{zk) - ZkZk-i)dui 

(143) 



k 



Let us introduce the multivalued function: 



rp 

Ui{p) := I dui (144) 



its discontinuities along the A and S cycles are: 

disc^^.Mj = , disc^^.Mj = 5ij (145) 
After an integration by parts, we have: 

4— = 2T?7i + 21^ + 2 tjTi^j - Res Ui{Vl,{zk)dzk - Zk+idzk - Zkdzk+i) 



+ y2 Res Ui{Vl{zk)dzk - Zk^idzk - Zkdzk-i) 

poo — 

k 

2Tr]i + 21^ + 2 ejTij - Res Ui{zk-idzk - z^dzk+i) 

3 k 

+ Res Ui{zk+idzk - Zkdzk-i) 

Poo — 

k 

2Tr]i + 2Ti + 2 > ejTjj — Res UiZ^idz^ + Res UiZj^dzf^j^i 
i 

+ Res UiZj\f^idzj^ — Res UiZ^dz^i 

Poo— Poo — 



= 2Trji + 21^ + 2 ^ ejTjj + Res UiZj^dzj^+i + Tui{poo-) - Res Wj^oc^^-i 

Poo+ Poo — 

J 

(146) 

using Riemann's bilinear identity and eq. p33|) . we find eq.( 



5.3 Derivatives with respect to T 

We have: 



^ = -ln77 (147) 
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indeed, from eq. (jll2|) . we have: 



df 



■zo dT 



20 



r 1 

)dZQ+ / ( — 



+1 



dT 



using eq. ()103|) and eq. (j^I?jl . we have: 



)dz^ — In 2;o — In z_\f 
(148) 



df 



P ^dzj^ dA 

+ -:-) - ln2;o - in z^f 



ZM A 



Zq 



In — — In 7 + In zj\fA — In 7 — In — In zj^ 



In 77 (149) 



We have: 



dF^^) A/- + 3 

:/i-T(— — +7V^ In T) 



(150) 



dT ' 2 
Indeed, from eq. fjllip . and using eq. (jl(J3j) . we have: 

dF du 
4— = 2ii + 2T^ + 2y"eirii + y"Resiyk{zk)-ZkZk+i)dS 
01 01 ^-^ ^-^ P00+ 

i k 

-y"Res{Vk{zk)- ZkZk-i)dS -AT{l+M\nT)-2MT{lhl) 



using eq . (11471) and eq. (jl33p . we find eq. (jl5(J|) . 

Remark: The derivative of the free energy wrt to T can be computed 
directely from eq.(ITI): 



d{F/T^ 



dT 



N 
2^2 



AT 



tr MMk) - MkMk-i ) (152) 



fc=0 



fc=i 



The matrix integral is defined for fixed fiUing fractions, i.e. for fixed ej/T. 
We can thus write: 



T 



,d{F/T^) 



dT 



f dF ^ I ^ ^ \ 

i * \ fc=0 k=l I 



Using the change of variable 5Mk = Mk, we get the loop equation: 

2T = ^ (tr MkV^iMk) - Mk{Mk-i + Mu+i)) 



(154) 
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which imphes: 

= i / tr X](y.(M,) - \M,Vi{M,))\ + {^^+l)T 

i \ k=0 J 

(155) 

i.e. 

F dF \ ^ 1 

2- - — - - ^ e,r, = Res (\4(Mfe) - -MfcV;'(Mfc))iyfcrfzfc + {M + l)T 

i fc=0 

(156) 

which is equivalent to eq. ()150|) . 

Note that eq. p54p , is nothing but the infinitesimal version of the rescaling 
Mfc — > aMfc. In particular one can choose a = a/T, and get directely from 
eq.©: 

F((?fc,,-, e„ T) = T2F((7fc,,T^/'-\ ~gk,jT^'^-\ e./T, 1) - ^^T^ InT 

(157) 

Taking the derivative of eq. (jl57j) with respect to T, gives again eq. (jl5(J|) . 



6 Other observables, leading order 
6.1 Resolvents, leading order 

So far, we have found that Wq{zq) and Wj^{zj^) obey algebraic equations. We 
have been able to determine the resolvents for the matrices at the extremities 
of the chain, but not for intermediate matrices. We are going to determine 
Wk{zk) for < A; < AT. 
We start from eq. ()132|) : 

W,{z) = --i- / ^^dz, (158) 
2i'kT Jc^ z- Zk 

where the integration contour is Ck defined in 14.31 Indeed, the residue is 
a contour integral around Poo+, and since equation eq. ()132|) was derived 
formally for large z (i.e. order by order in the large z expansion), we have 
to assume that the contour of integration encloses all the solutions of 
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Zk{p) = z. Moreover, the residue is the sum of poles at the p+j^fc(zfc)'s and at 

Poo+'- 



V"{zk-i)dzk-i - dzk-2 



Ck Vk-li^k-l) - Zk-2 

'*= 1 /" Viz )z 



Zk-i{p+j,kizk)) ^ '-^^ (159) 



9k-l,dk-i+l 



Therefore: 



9k-i,dk-i+i 9k+i,dk+i+i 

(160) 

Notice that we have: 

rk 

£ Zk.^{p,A^)) = SkKiz) + J^^±l^ - ^'-^''^-^ (161) 

9k+i4k+i+i 9k-i,dk-i+i 

6.2 2-loop functions, leading order 

From eq.dHZI) and eq. fTTpl We find {I < k): 



Tr — Tr — ) dzkdzi = -J^Y] B{p+i,k{zk),p+j,i{zi)) 

Zk- Mk zi~ Ml , , 



(162) 

6.3 2-point one loop functions 
6.3.1 Wk,k+i 

Define the following polynomial of two variables: 

Qk,k+i{,Zk, Zk+i) = Res dzQ... dzk-idzk+2 ■ ■ ■ dzj^ ' " " " ' 
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— E{Z_Q, . . . ,Z_i^_i, Zk, Zk+l,Z_i^+2^ ■ ■ ■ i^Af) 

fe-1 Af rk 

j=0 j=k+l 0^j=-Sk 

k M rk+i 

= (-1^+^ n^ivi n n (^^c-^fefe^+ii^fe+i))) 

j=0 j=k+2 Q^j = -sk + i 

(163) 

where the Z_j{zk, Zk+i) are defined by: 

Zj^^{zk, Zk+l) = Vi'{Z_j{zk, Zk+i)) - Z_j_^{zk, Zk+i) j > k + l (164) 
Z_j^-^{zk, Zk+l) = Vi'{Z_j{zk, Zk+l)) - ^j+i(2;fc, Zk+l) j < k 



We conjecture (proved in appendix B for A; = 0): 

1 — Wk,k+i{zk, Zk+l) = 



Qfc,fe+i(^fc, Zk+l )/n =0 fi'i.dj+i 1 lj=fc+2 9j,dj+i 



U%ii^k+l{p+j,k{Zk)) - Zk+l) Ilf=liZk{P-j,k+liZk+l)) - Zk) 

Tk ni=Y(^fc ~ Zk{p+j^k+i{zk+i))) 



(165) 

Remark: 



9k4k+i n;Li(jfe+i(p+.,fe(^fe))-^fe+i) 

Sfe+1 11^=1 ~ ^fe+l(P-j,fe(^fc))) 

Y[7=l\^k{P-,,k+l{Zk+l)) - Zk) 



Qk,k+i{zk, Zk+l) = Qk-i,k{Vl{zk) - Zk+l, Zk) (166) 

and we may conjecture that the spectral curves of the differential systems 
defined in the appendix of are: 

det (zk+ilrk+sk - Vkizk)) = (-1)- ^_^^w(^fe^^fe+i)^^ (ig7) 

1 [j=o 9j,dj+l 1 Lj=k+1 9j,dj+l 

so that the property eq. ()166p would be nothing but the duality discovered 
in 0. 
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6.3.2 The function U in the large hmit 

Since the function U appears in the RHS of the master loop equation, it is 
important to be able to compute it in the large N limit. 
Define (notice that Ui = U): 

k-l 

Uk{zo, Zk, . . . , z^f) := Pol W{zo, . . . ,zj^) fk,Af{zk, ■ ■ ■ ,zj^) Y\ Zj (168) 

Zl,...,Zj^ 

We shall prove that: 
Ukizo, Zk,..., zj^) = Hk,M{zk, zj^)Wo{zo) 

l=k J^k,l-l{Zk, Zl-i) 

I Y^AT-l E{Zo,...,Zi_^,zi,...,zj^)~E(Zo,...,Zi,zi+i,...,zj^) ^ / x 

+ 2^i=Q jpzi Hk,i-i{Zk,...,Zi-i) 

(169) 

where 

Hk,i{zki ■ ■ ■ , Zi) := Pol zo,...,zj^ 

j=0 ^ ^ j<k j>l 

i.e. Hk^i is a polynomial in Zk, ■ ■ ■ , zi, and satisfies: 

Hk+i,k = 1 , Hk^i = if A; > / + 1 , Hk^i = ^ — ^-Hk+i^i — Hk+2,1 

Zk — ^k 

(171) 

In particular for = 1 eq. (jl69|l reduces to: 
V {zq, ...,zm) = V^{zm) - zat-i 

S-^J\f-l E{Zo,...,Zk,Zk+i,--;ZM)-E{Zo,...,Zk+i,Zk+2,---,zj^) tt ( ^ ^ \ 

+ Z.fc=o tLi^k[zu...,Zk) 

= Hi^x{zu . . . , zj^)Wo{zo) 

EA^-l P{Zo,...,Zk,zk+l,■■■,z^r)-P{Zo,...,Z^:+l,zk+2,■■■,z^r) tt r ^ ^ \ 

k=0 h,,k{Zi, . . . , Zk) 

(172) 

Proof of eq.dinnD: 

Let us define: 

k-l 



Ak{zo, Zk,..., zm) := Pol TT Zj fk+i,Ar{zk+u . . . , zm)W{zo, . . . , zj^) 

Zi, 1 .2; J. I 1 z\r 

(173) 



Zl,...,Zk-l,Zk + l,---,Zj>/ 
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Sk{zo-, Zk-, ■ ■ ■ ■, zj^) :— Pol fo,k-i{zo, . . . , zi:_i) 

ZO,...,Zk-l,Zk + l,...,Zj^/ 

fk+iMi^k+1, • • • , zj^)W{zo, ...,zj^) (174) 



Tk{zo, Zk,..., z^f) ■= Pol fo,k~i{zQ, Zk^i)V^{zk) 

zo,...,z^ 

fk+i,M{.Zk+i, zj^)W{zo, ...,zj^) (175) 

From eg . (1131) we have: 

Uk = Pol V^Ak - Uk+2 (176) 

and from eq. ()29|) . we have: 

{zk — Zk)Ak = Uk+i — Sk{Zo, . . . , Zk-i, Zk, . . . , zj\f) (177) 
That implies: 



_ V^{zk)-VaZk) 

'Jk — 7^ t/fc+l — Uk+2 

Zk — ^k 

TkjZQ, ■ ■ ■ , Zk-i, Zk, ■ ■ ■ , zj^) — Tfc(Zo, ■ ■ ■ , Zk, Zk+i, • • • , ^^^ -^jQ-j 

Zk — Zk 

It is easy to prove (by recursion), that /^at — fo,k-iVl{zk) fk+i,Af is hnear 
in Zk, in other words we can write: 

hM = k,k-iVl{zk)fk+i,M + ZkBk (179) 

where is independent of z^, therefore: 

Tk{zQ, . . . , Zk, . . . , zj^) - P{zo, ...,Zk,...,zj^)= independent of z^ (180) 

and thus eq. ()178|) can be rewritten: 

VlXzkl^ViiZkl 

Uk — Uk+1 — Uk+2 

Zk — ^k 

P{Zq, . . . , Zk-i, Zk, ■ ■ ■ , zj\f) — P{Zq, . . . , Zk, Zk+i, ■ ■ ■ , ^A/")^ -j^g-j^-^ 

Zk — Zk 

From which, the initial conditions Uj^+i = Wq, = (easily derived) 

imply eq. (jl69|) . 
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6.3.3 The extremities correlator Wq^^ 
From eq. ()172p and eq.dH}, we derive: 

W,M.«, ^U) = E ^'^'-■■y^^- (182) 

where: 

Zk ■■= Zk{p+ifl{zo)) , Zk := Zk{p-i,Af{zj^)) (183) 



7 Subleading epxansion 

The aim of this section is to generalize the calculation of [llUlllj . and compute 
the next to leading 1/iV^ term in the topological expansion. In this purpose, 
we expand the various observables in a l/N"^ power series: 

Zkizo) = Z^^\zo) + ^X^'o) + ■ . . (184) 

P{zo, ...,z^) = P(°)(^0, + -^P^'\zo, ...,zj^) + ... (185) 

and so on. Then we expand eq.()28p to order l/N"^: 

X: Z^^\z,)Ek - P^'\zo, Zr\ . . . , Z;?)) = f/(^.o, Zr\ Z^J); z^) (186) 

k=l 



where: 

dE(^\zo,...,z^f) 



Zi=Z 



(0) 



(187) 



dzk 

From the definition of the Z^s eq. fjlH) . we easily find: 

zl^\zo) = Hi,k-iZ['\zo) (188) 
where (the Hk/s have been defined in eq. (|17(J|) ): 

:=/fi,,(zf ,...,zf) (189) 

We have: 

= , H,,o = 1 , Hi,k = V;:{zf^)Hi,k-i - Hi,k-2 (190) 
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I.e. 



Therefore eq. 



reads: 



d''fi,k{zi, . . . ,Zk) 



^z^ ... dzu 



(0) 



(191) 



Y^^=l Hl,k-l Ek 



(192) 



Similarly to what was done in ^HIEI; we determine the polynomial P'^^\ 
by the condition that z'^\zq) has singularities only at the endpoints eo,i 
in the zo^sheet. That condition is sufficient to determine all the unknown 
coefficients of P^^\ Now, let us compute U^q"* ■ 



7.1 Computation of U-q{zq, . . . , zj^f, x) 
Notice that there is no explicit depedence on Vq in eq.()26|). therefore: 

dU{zQ, ...,zm) 



U-o{zo, . . . , zj^; z) 



dV.iz) 



(193) 



For /c > 1, write Zk = zj^^ + C^, and Taylor expand eq. ()172|) in the C's: 
U {zo, ...,zx) = Vj^{zj^) - zx^i 

Af ^ AT 

- + 7^(kEk,k + ^ CiEk,i) 



k=l 



i=k+l 



fc-1 



[Hl^k~l + CiHl,k-l;i) 



(194) 



1=1 



where Ek was defined in eq. (ll87p . and Ek^i is: 

d'E('^\zo,...,ZM) 



E, 



k,l ■- 



dzkdzi 



z^=Z. 



(0) 



(195) 



H\^k was defined in eq. p89p and: 

dHiJzi, ...,Zk) 



Hl,k;i 



dzi 
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(0) 



(196) 



From eq. (|17Up and eq. (|189|) we derive: 



(197) 



Then take the d/dVo{z) derivative, using d(k/dVo = —dZk/dVo, and take 
the C hmit: 



Ufi^ i^o, zf^ , . . . , Z^^ ; z) — ^ 



dE, dzf E,,, ^ dzf 



^ ( dHik-i dZ^^^ _ 



k=l 



(198) 



From eq. (fTT|) . one easily derives: 



dzit dZ 



dzf\ - dZ\ 



dV^{z) dV^{z] 
and from eq. ()189|l 



v;"(^ro - 



(0) 



dV^{z) ''^dV^{z) 



(199) 



(200) 



7.2 Computation of dEk/dVQ{z) 
Let us define Zj{zo,Zi): 

Zq{zq,Zi) := Zq , Zi(zo,Zi) := Zi , Zj^i(zo, Zi) = Vj{Zj(zQ, Zi))—Zj_i(zo, Zi) 
We have: 



(201) 



dZj{zo,Zi) 



dzi 



Hi 4^1 , 



d'^Zk{zo,zi) 



k-l 



Consider the polynomials defined in eq . (11631) : 

Qo,iizo,Zi) := E^'^\Zo{zo,Zi),...,Zx{zo,Zi)) 



(202) 
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so 

C(^i - Zf\z,)) n(^i - zi(p-,,o(^o)))) (203) 



(where C = — Iljli fi'j.dj+i)' take its derivative (using eg. ()2()2|l ): 



dQoAzo, zi) 



dzi 

That imphes: 



J2EjH,,_, = Cn(^r - ^i(P-.,o(^o))) (204) 



so 



5^i(p_,-o(^o)) 



Zf^ - 2;i(p-/,o(^o)) 



- \dVo{z) dVoiz) 
Plugging that into eq. ()198p . we get: 

^° ^ dzi'^ azi(p_,,o(^o)) 



(205) 



(0). 



J-1 



E 



so 



- zi{p-i^o{zo)) 



dV^{z) 



f dz['^ g^i(p_,,o(^o)) 



5K)(^) 



Zf^ - zi{p^ifi{zo)) 



zo. 



dZ^ ^ V ( Ek k 2 \ 



so 



-E 



92;i(p_j-o(2;o) 



dV,{z)) 



zo 



(E^Ii ^1 
Zi'^ - Zi{p_ifl{zQ)) 
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1 ^ ^ — — 

-2 9l^5ZE^^.^-^i''^-i^i.^-i 



1 dZ\ 



k=i j=i 

(0) AT fc-1 



2 9K)(^; 



(206) 



fc=i i=i 

Take the second derivative of eq. (|2U3|) : 



5^Qo,l(-2o,-2l) 



9zf 



so 



zi=Z\ 



(0) 



AT AT 



fc=l J=l fc=l 1 

plugging eg . (12071) and eq. (j202p into eg . (12061) . we get: 



(207) 

Zl=Zl 



(0) 



7(0). .^ 



1 



dzi{p-jfl{zo)) 



so 

lr^--i(p-/,o(^o)) ^^o(.) 



^0 



(208) 



which can be compared to [TT 



7.3 Next to leading order 

Using eq. ()208p into eq. ()192p . and using eq.(|H7j). we get the correction 
for the resolvent: 



z['\zo) 



1 



dzi {p_jfi{zo)) 



so 



lr^r--i(p-.,o(^o)) ^^o(.o) 

5(P-«,o(^o),Pi,o(2:o)) 
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^ 2;i(pi,o(2o)) - Zi{p-i^q{zq)) dzo{pi,o{zo))dzQ{p_i^oizo)) 
(0) 7(o)x 

f ^^■■■^y> (209) 



X]fc=i Hi k-i Ek 
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and we remind that P^^^ is completely determined by the condition that 
z[^\zq) has singularities only at the endpoints eo,i in the ZQ-sheet, and has 
vanishing i3-cylce integrals. 

From there, it should be possible to extend the calculation of [IHIIII] to 
the chain of matrices, and compute the term in the free energy for the 
chain of matrices. This will be left for a later work. 



8 Examples 

8.1 Example: one-cut asumption (genus zero) 

Let us assume that the genus is 5^ = 0. 

It was already discussed in [10] that for multimatrix models, the so called 
one-cut asumption should be replaced by a genus zero asumption. Indeed, 
in that case, the number of zo^endpoints, i.e. the number of zeroes of dzQ 
is equal to the number of sheets, therefore there is exactly one cut in the 
physical sheet. S is in one to one correspondance with the complex plane 
C, and it is well known P^ ITTj that the parametrization eq.()81|) is rational, 
and with s = it can be written: 

P-Pco+ ' 



(210) 



Zk{s) has two poles: one pole of degree at s = 00, and one pole of degree 
Sfc at s = 0. This parametrization is identical to the one found in ^3] by the 
biorthogonal polynomial's method. 
The set of equations: 



Vl{zk{s)) = Zk+l{s) + Zk-l{s) 

V^izois)) - z,is) 



T 



Ts 



(211) 



"0,1 = "AT -1 =7 = 7 



is sufficient to determine all the a^^j (one is free to impose ao,i = (yj\f-i 
because s can be changed into any constant times s). 
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The Abelian differential of the third kind is: 

ds 



dS = — 



The Bergmann kernel is: 



ds ds' 
(s - s'Y 



(212) 



(213) 



8.2 Example: gaussian case 

Consider all potentials quadratic (i.e. dk = 1) 

T/ / \ 9k 2 

yk{zk) ■■= 



(214) 



A direct computation of the matrix integral eq.Q gives the free energy: 



2 2 



(215) 



Define: 



Dki := det 



f9k 1 
1 ■•• 

V 



■■. 1 

1 9iJ 



Dk,k-i := 1 , Dk,k-2 := (216) 



The function fk,i of eq. (fT^ is: 

fk,l = Zk . . . Zi Dk,l 
thus the polynomial P is a constant: 

P{zo, ...,Z^f) = 

and the polynomial E is: 

E{zo, ...,z^f) = {gozo - zi){gj^fzx - ^at-i) - TP 
The leading order loop equations are thus: 

Zk+l + Zk-l = gkZk , ^-l^AT+l = TP 



(217) 



(218) 



(219) 



(220) 



47 



Note that eq. ljH^ implies that the genus is necessarily g = 0, and thus 
there is a rational parametrization of the form eq. (j21U|) . with = Sk = 1, 
namely: 

^fc = J TT— (^o,fc-is + Dk+i,^s-^) (221) 

We find: 



7 = 7=^7^ (222) 
and: 

/i = r + rinDox-^lnr = T-rin77 (223) 
and it is easy to check that eq. piip coincides with eq. ()215j) . 

8.3 Loop equations of the 2 matrix model 

Let us take T = 1. The 2-matrix model corresponds to A/" = 1 we have: 

/o,i(^o, ^i) = K)'(^o)K(^i) - ^0^1 , /i,i(^i) = K(^i) (224) 

i.e. 

and the master loop equation reads: 

(Fo'(^o) - Z,){y[{Z^) - Zo) - P(Zo, Zi) = ^f/(^, ^i; Z,) (227) 
with Zi = Vq{zo) - Wq{zq). We recover the equation of jl2H13[ITnj . 

8.4 The 3 matrix model 

Let us take T = 1. The 3-matrix model corresponds to A/" = 2, i.e.: 

/0,2(^0, Zi, Z2) = V^iZo)Vl{z^)V^iz2) - ZoZiV^{z2) - Z2ZiVl{z^) (228) 



48 



I.e. 



fiM, Z2) = Vl{z,)V^{z2) - z,Z2 (229) 



1 / V^jzo) - V^jMo) V;{z^) - VljM,) Vi{z2) - ^^M^) \ 



2A^ \ ^ zq-Mq zi-Mi Z2-M2 I 

1 / V!,{Z2) - Vi{M2) Vi{z^) - VliM,) V^izo) - V^{M^) 



(230) 



2A^ \ Z2- M2 zi - Ml zq- Mo 

1 / Ki^o) - ypiMo) \ 1 / V^iz2) - vm2) \ 
N\ zo~Mo / N\ Z2-M2 I 



Uizo.z..Z2,z) ^ 1/tr ' Vi{z.)-Vi{M,)Viiz2)-ViiM2)^^ 1 



2 \ ;zo - Mo zi- Ml Z2-M2 z- Mo 

I 1 / ^2 (^2) - ^^(M^) V;{zi) - VljMi) 1 1 



2 \ Z2-M2 zi- Ml zo- Mo z- Mo 

/I 1 
- ( tr — tr — 

\ zo-Mo z-Mo 

(231) 

and the master loop equation is: 
(Vo'(Zo) - Zi)(l^2'(^2) - ^1) - P{Zo, Zi, Z2) = ^U{Zo, Zi, Z2; Zo) (232) 

with Zo = zo, Zi{zo) = Voizo) - Wo{zo) and ^2(2:0) = V{{Zi) - Zo. 

Large one loop functions 

We define: 

Zi{z2) = Vi{z2) - W2{z2) , Zo{z2) = ^(^1(^2)) - ^2 (233) 

1^6 1(^0, zi) = Pol V^{zo)Wo,i{zo, zi) (234) 

2(^1, ^2) = Pol Vi{z2)Wi.2{zi, Z2) (235) 

1 2(^0, zi, Z2) = Pol K,'(^o)V^2(^2)W^(^o, zi, Z2) (236) 

' ' 20,22 

eq.dlHS) reads: 

M.o.(.o, = —V f - ^(^^ I (237) 

Zi — Zi \ 2:0 — Zo Z2 — ^2 
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eq. (|169|) reads: 

U2{zo, Z2) = Pol V^{Z2)W0,2{Z0, Z2) = ^^^'^ ^'^^'^ (238) 

Z2 Z2 — Zj2 

eg . (11721) reads: 

'U{zq^Z\^Z2) = Pi^2(yZo, Zi, Z2) 

-.ri, X , E{zq,Z\,Z2) — E{zq,Zx,Z2) 
— y2[Z2) - Zi -\ z - Z 

, VI{z,)-VI{Z,)E{zoX,Z2) 

+ ; — ^ ; — ^ — l^-^y) 

Zi — Z-i Z2 — ^2 



eq. (jH2|) implies: 

o ^1 (^1) - '^1 (^1) u/ / N A,2(^o, ^1, ^2) - ^1,2(^0, ^1, 2:2) 

— yyo,2[zo,z2) 

Zi — Zi ^ Zi — Zj\ 

(y(-Zo-Z2)( E{zo,Zi,Z2) E{zo,Zi,Z2) 



{Z,-Zi) \{zi ~ Zi){zo - Zo) (z, - Z,){z2 - Z2) ^ 

+ 1 E{zo,z„Z2)^ ^240) 

{zi - Zi){zi - Zi) 

We have the relationships: 

{z, - Zi)Po,o(^o, zi, Z2) = ^^'^'^'^''^ - W(,^,^^{zo, zu Z2) (241) 

— Zq 



ir. f^^.rrl X ^(^0,^1,^2) ^(^0,^1,^2) 

[Zi- Zi)W[Zq,Zi,Z2) ^ 



{Z^ - Zi){zo - Zo) (zi - Zi) {Z2 - Z2) 

{Zi - Zi)Wo 12(^0, 2:1, 2:2) 



{zi - Zi){zi - Zi) 
But so far, we have not been able to compute PVq We conjecture: 



(242) 



(2^0 + Z2- VI{zi))Wq 1 2(^0, zi, Z2) = E{zo, zi, Z2) 

-(K,'(^o) - ^1 - 1^6,1(^0, ^i))(V^2(^2) - ^1 - 1^1,2(^1,^2)) (243) 

indeed, both sides are polynomials in zq and Z2 with the same degree and 
same large zq^ z\ and Z2 behaviours. 
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9 Conclusion 



We have written the loop equations in a very exphcit way. To leading order 
the loop equations become algebraic, and using that fact, we have derived 
many observables, like the free energy and many correlators. 

We have computed some mixed correlators and we have given conjecture 
for others. So far, the complete one-loop function has not been computed, 
and we don't even have a conjecture (except in the 3- matrix model). The 
mixed correlators are important in the study of boundary operators [20] • 

It would also be interesting to generalize ^U] to the chain of matrices, 
and find the next to leading large N corrections to the free energy. 

Aknowledgements The author would like to thank V. Kazakov for encour- 
aging him to translate in english (and add new results to) old notes written in 
french, and to thank M. Bertola for fruitfull discussions and careful reading. 
When this work was in progress, M. Bertola found independent proofs for 
some of the formulae in section 5, which will be published later. 



Appendix A Some usefull formulas with residues 
and polynomial parts 

We make an abundant use of the following formulas: 

p^m^y'(')-nz) ^ i::(£M£ . _r(z) (2) 

z z — Z z — Z z — Z 

V'{z) ^ _J_ ( V'{z)-V'{M) _ V'{Z)~V'{M) \ 
z {z-Z){z-M) z-Z\ z-M Z-M ) 

V'{z)dz ^ V'{Z) - V'{M) 
{z~Z){z-M) Z~M 

V'jM) ^ V'{z) _ V'jz) - V'jM) ^ V'jz) _ V'jz) 

z-M z-M z-M z-M z z-M ^ ^ 

M z ^ z ^ ^ z 

1 = Po (6) 

^ i\/r ^ ^ i\ /T \ ' 



Z-M z-M z-M z z-M 
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Appendix B Determination of Wo,i 

Define the following polynomial of zq: 

W(,,{zo,zi) = Pol V^{zo)Wo.i{zo,zi) (1) 

eq. ()H2j) implies the following identity (to large leading order): 

{z, - Zi)(l - 1^0,1(^0,^1)) = zi- V^{zo) + W^o,i(-2o,^i) (2) 
which can also be written: 

Z-i VJ.(zo)-WnJzo,zi) 
(1 - ^)(1 - WoA^o, z,)) = 1 - ' (3) 

Zi Zi 

Notice that the RHS is a polynomial in zq. 
take the log: 

^1 y^izo) -Wn,{zo,Zi) 
In (1 - ^) = - In (1 - WoA^o, z,)) + In (1 - ^ ) (4) 

Zl Z\ 

take the fractionary part (Frac/ := / — Pol/), we find: 

-ln(l-W^o,i(^o,^i)) = Frac (5) 

20 \ 2:1 y 

i.e. 

W^o,i(^o,;^i) = l-exp(^-Frac (^In (1 - ^^^)^ ^ (6) 
in other words: 



^ ? rf^i(p) ^^.^ ^o(p). 
2i7r 7ci ^i(p) - ^1 Zq ' 

Frac f In (1 ^^^"^^^ 



21 V Zq 



In (1-^^)- Pol fin (1-^5^)^ (7) 

ZO ^1 V Zq ) 



where the second equality is obtained by integration by parts. 
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In other words we have that: 

(1 - Wo^i{zo, zi)){zi - ^1(2:0)) = Polynomial in 2:0 



(8) 



(1 - 1^0,1(^0, 2:1)) JJ^ (2:0 - Zo(p_j,i(^i))) = Polynomial in 2:1 (9) 
that implies that: 



1 - Wo^i{zo,Zi) 



Qo,i{zo, Zi) 



{zi - Zi{zo)) Yl7=ii^o - Zo{p-j,i{zi))) 



(10) 



where Qo^i{zo, Zi) is some polynomial in both variables, of degree Vi + Si, ro + 
So, which vanishes each time there exists p such that zq — Zo{p) and zi — 
Zi{p). That implies that: 



ri 



Qo,i{zo, zi) oc Ylizo - ^o(p+j,i(^i))) n^^o ~ ^o(P-j,i(^i))) 



ro 



so 



(X Ylizi - Zi(p+j-o(^o))) n(^i ~ ^i(P-i,o(^o))) 

3=1 j=l 

(11) 

i.e. (5o,i must be equal to (same degree, same zeroes): 

E{zo, zi, 2^2, ... , z^f) 



I E 
(9o,i(^o, zi)(x (b dz2... dz^f 

J 1 i,=i 



{zj+^ + Zj^i-Vj{zj)) 



(12) 



We thus have: 



Wo,i (^0,2:1) = 1 



<5o,i(^o, Zi) 



{zi - Zi{zq)) nyLi(^o - ^o(P-j,l(^l))) 



(13) 



I.e. 



W^o,i(^o,2;i) = 1 +5^0,^0+1 



]\%i{zq - Zo(p+j-i(zi))) 
{zi - Zi(p+i,o(2;o))) 



(14) 
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